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Abstract. Association schemes were originally introduced by Bose and his 
co-workers in the design of statistical experiments. Since that point of incep- 
tion, the concept has proved useful in the study of group actions, in algebraic 
graph theory, in algebraic coding theory, and in areas as far afield as knot 
theory and numerical integration. This branch of the theory, viewed in this 
collection of surveys as the "commutative case," has seen significant activity 
in the last few decades. The goal of the present survey is to discuss the most 
important new developments in several directions, including Gelfand pairs, co- 
metric association schemes, Delsarte Theory, spin models and the semidcfinitc 
programming technique. The narrative follows a thread through this list of 
topics, this being the contrast between combinatorial symmetry and group- 
theoretic symmetry, culminating in Schrijver's SDP bound for binary codes 
(based on group actions) and its connection to the Tcrwilligcr algebra (based 
on combinatorial symmetry). We propose this new role of the Terwilliger al- 
gebra in Delsarte Theory as a central topic for future work. 



1. Introduction 

The concept of (symmetric) association schemes was first given in the design of 
experiments |27l I28| . It can also be viewed as a purely combinatorial generalization 
of the concept of finite transitive permutation groupsQ The Bose-Mesner algebra, 
which is a fundamental tool in the theory, was introduced in [26j . The monumental 
thesis of P. Delsarte proclaimed the importance of commutative association 
schemes as a unifying framework for coding theory and design theory. In [52] . 
the method of linear programming was successfully combined with the duality of 
the (commutative) Bosc-Mesner algebra, and it has been serving as one of the 
underlying principles in the theory of commutative association schemes. 

The theory continues to grow rapidly, enlarging its diverse connections and ap- 
plications to other branches of mathematieso However, a number of important 
topics are either just glanced upon or not treated at all in this paper. We are 
limited not only by length, but also by our own limited expertise. Notable omitted 
areas include: distance-regular graphs^ the Terwilliger algebraQ and the role of 



2000 Mathematics Subject Classification. 05E30. 

Key words and phrases. Association scheme; Delsarte Theory. 

1 The theory of association schemes is famously said to be a "group theory without groups" 
[20]. 

2 Bannai 1111 p. 108] also states that "it seems that commutative case forms a fairly closed 
universe, similar to the compact symmetric homogeneous spaces." 

3 Indeed, a separate, comprehensive update on distance-regular graphs by different authors is 
reportedly in preparation and was originally planned for this collection of surveys. 

4 While we cannot begin to cover this important topic in the present paper, we shall encounter 
the Tcrwilligcr algebra, particularly in our vision for Delsarte Theory and our treatment of semi- 
definite programming. 
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association schemes in designed experiments. In fact, if we had tried to cover all of 
these topics in full detail, then this survey article would have been a book rather 
than a paper! However, fortunately there are a number of excellent articles/books 
on the theory of (commutative) association schemes and Delsarte Theory, e.g., 
[52 [201 EQl E3 [311 HI HO- The following are a few of the books which include 
accounts on commutative association schemes: |110[ l35 | Ff3 l 1104] . Thus, naturally 
guided in part by own current research interests, we shall in this paper focus on 
recent progress in the theory that has not been treated (in detail) in previous lit- 
erature. Thus our hope is to contribute an addendum to the important references 
listed above; in no way do we intend to supplant any of them or summarize their 
content. 

We now give a summary of the present paper, together with some additional 
comments. Sj2] reviews the basic theory of commutative association schemes, the 
goal of the exposition being to provide just enough background to support the later 
sections. We refer the reader to [20l [30] for a more comprehensive and detailed 
account of the theory. In Sj3]we consider commutative association schemes obtained 
from transitive group actions. The Bosc Mcsncr algebra of such an association 
scheme coincides with the corresponding ccntralizcr algebra, and its eigenmatriecs 
are equivalent to the zonal spherical functions. We briefly introduce a beautiful 
phenomenon which can be observed in many families of commutative association 
schemes related to classical groups over finite fields, and discuss its recent applica- 
tion to the construction of Ramanujan graphs which are of great interest in com- 
puter science. The topics of S|H are metric (or P-polynomial) association schemes 
and cometric (or Q-polynomial) association schemeso These concepts were intro- 
duced by Delsarte |52j . and certain systems of orthogonal polynomials naturally 
come into play in the theory. The association schemes which are both metric 
and cometric may be viewed as finite analogues of rank one symmetric spaces |20[ 
pp. 311-312], and many researchers have been working on the classification of such 
association schemes. Special attention will be paid, however, to the class of comet- 
ric (but not necessarily metric) association schemes, which has become an active 
area of research in its own right over the last decade. 

§ [U discuss Delsarte's theory and related topics. We especially recommend 
the survey articles [59l [37] for the (major) progress up to 1998. We shall see that 
codes and designs are dealt with in a unified manner within the framework of com- 
mutative association schemes. While the minimum distance and dual distance (or 
maximum strength) are two important parameters of a code or design in Delsarte 
Theory, a similar theory exists — with an entirely different class of examples — 
for two new parameters, namely width and dual width, which were introduced by 
Brouwer et al. in [31j . These parameters will also be briefly reviewed in $5] In ^J6j we 
first recall standard facts on translation association schemes and their group codes. 
The duality between codes and designs observed in $5] becomes more than formal in 
this case. We also discuss dualities of association schemes in connection with spin 
model theory (see e.g., [90]). The subject of fj7] is the famous linear programming 
bound of Delsarte. Delsarte himself explored the specialization of this technique 
to metric/cometric association schemes in |52j . and there are a number of excellent 
treatments on the impact of this technique in coding theory (e.g., [1101 IT03j). Our 



The concept of metric association scheme is essentially the same as that of a distance-regular 
graph, which is in turn a combinatorial version of the concept of a distance-transitive graph. 
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update takes as its point of departure the striking 2001 result of A. Samorodnit- 
sky |141j which says that something beyond the linear programming bound will 
be required to resolve the most fundamental problem in algebraic coding theory, 
namely the determination of the asymptotically optimal rate of a binary block code 
for a communication channel with given bit-error probability. This serves as strong 
motivation for what is to come in §101 But we also aim to extend the techniques 
beyond the metric/cometric cases in support of the applications in <j8] Dclsarte's 
theory has been quite successful for codes/designs in metric/cometric association 
schemes, but the purpose of $S]is to introduce far broader applications of his theory. 
Working mainly on the "eigenspace side," we demonstrate the great value in extend- 
ing the theory beyond the class of cometric schemes by simply listing a variety of 
unusual settings where Delsarte's theory applies. That is, we shall characterize var- 
ious combinatorial objects as codes or designs in certain association schemes. Here, 
the linear ordering of cigenspaces fundamental to a cometric association scheme is 
replaced by a partial order and certain well-known posets play a key role in the 
study of designs, and in finding solutions to the linear programming bound. For 
example, (t, m, s)-nets |137( HIT] — which provide quasi-Monte Carlo methods for 
numerical integration, simulation and optimization — are closely related to the Del- 
sarte designs in the "ordered Hamming scheme," a family of association schemes 
generalizing the ordinary Hamming scheme. Here the eigenspaces are indexed by 
the members of a downset (or "lower ideal") in Young's lattice. 

In fj9l we briefly investigate the Terwilliger (or subconstituent) algebra of an 
association scheme |163[ I164| 1165] . This noncommutative matrix algebra contains 
much more structural information about the association scheme than the (commu- 
tative) Bose-Mesner algebra. The Terwilliger algebra has proven to be a powerful 
tool in the study of metric & cometric association schemes. Despite the importance 
of this connection, our focus in this section is instead on the use of this algebra in 
the analysis of codes and designs, mirroring the use of the Bose-Mesner algebra 
in Dclsarte Theory. This is still quite a new approach, but we have included the 
account here to propose it as a possible research direction. CTUl is devoted to a dis- 
cussion of the semidefinite •programming bound introduced recently by A. Schrijvcr 
|143| . This new bound was first established for binary codes and constant weight 
codes in |143] using the Terwilliger algebra, and then for nonbinary codes in [70j . 
The semidefinite programming bound is always at least as good as Delsarte's linear 
programming bound, and numerical computations show that there are many cases 
where it in fact improves upon known upper bounds. For simplicity, our exposition 
is restricted to binary codes. A survey paper [173] also contains an account on this 
bound for binary codes based on the results in [62], but we shall particularly em- 
phasize a viewpoint centered on the Terwilliger algebra. The current formulation 
of the semidefinite programming bound relies heavily on certain group actions, so 
that strictly speaking it does not belong to the "association scheme theory" yet. 
However, experience shows that group actions can be supplanted with assumptions 
of combinatorial regularity and our treatment here reflects our hope that, for a wide 
class of schemes, Delsarte Theory will be reconstructed entirely in the near future 
based on the Terwilliger algebra and the semidefinite programming method. 
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2. Commutative association schemes 

Let X be a finite set and C XxX the set of complex matrices with rows and 
columns indexed by A. Let 1Z = {Ro, Ri, ■ ■ ■ , Rn} be a set of non-empty subsets 
of X x X. For each i, let Aj g C XxX be the adjacency matrix of the graph (X, Ri) 
(directed, in general). The pair (X,7Z) is an alsocaUon schem® mtl n classes if 

(AS1) Aq = I, the identity matrix; 

(AS2) J27=o Ai = J, the all ones matrix; 

(AS3) Aj £ {A h A!,..., A n } for O^i^n; 

(AS4) AiAj is a linear combination of Aq, Ai, . . . , A„ for ^ i, j ^ n. 
By (AS1) and (AS4) the vector space A spanned by the A is an algebra; this is the 
Bose-Mesner (or adjacency) algebra of (X, TZ). The Aj are linearly independent 
by (AS2) and thus form a basis of A. We say that (X, TZ) is commutative if A is 
commutative, and that (X, TZ) is symmetric if the Aj are symmetric matrices. A 
symmetric association scheme is commutative. Below are listed a few examples of 
(symmetric) association schemes: 

Example 2.1. The Johnson scheme J(v,n) (v 2n): A is the set of all n-element 
subsets of a fixed set f2 with u points, and (x, y) £ Ri if \x fl y| = n — i. 

Example 2.2. The Hamming scheme H(n,q) (q ^ 2): A is the set of all words of 
length n over an alphabet Q with q symbols, and (x, y) € i?; if x and y differ in 
exactly i coordinate positions. 

Let C x be the set of complex column vectors with coordinates indexed by X, 
and observe that C XxX acts on from the left. For each x G X let x be the 
vector in C x with a 1 in coordinate x and elsewhere. We endow C x with the 
standard Hermitian form (, ) so that the x form an orthonormal basis for C^. 

For the rest of this section, let (A, 1Z) be a commutative association scheme with 
adjacency matrices Ao,Ai,...,A„ and Bose-Mesner algebra A. By (AS3) A is 
closed under conjugate transposition. Since A is commutative, it follows that there 
is a unitary matrix U £ C XxX such that U^ 1 AU consists of diagonal matrices only. 
In other words, C x is decomposed as an orthogonal direct sum of n + 1 maximal 
common eigenspacetQ of A: 

C x = V ±Vi±--- ±V n . (1) 

For each i, let Ei £ C XxX be the orthogonal projection onto Vi. Then the Ei 
form a basis of the primitive idempotents of A, i.e, EiEj = SijEi, y,"_n Ej = I. 
Note that |A| _1 J is an idempotent in A with rank one, hence must be primitive; 
we shall always set Eq = |A| _1 J. It also follows from (AS2) that A is closed 
under entrywise (Hadamard or Schur) multiplication, denoted o. The A form a 
basis of the primitive idempotents of A with respect to this multiplication, i.e., 
Aj o Aj = SijAi, X^i=0 = J- 

The intersection numbers p\^ and the Krein parameters g*j- (0 ^ ^5 n) of 

(A, 72.) arc defined by the equations 

n n 

M-Ep^> ^°^- = rx|E4^- ( 2 ) 

k=0 ' ' fc=0 



^The notion coincides with that of homogeneous coherent configuration; see 1641 . 
^They are sometimes called the strata [8]. 



COMMUTATIVE ASSOCIATION SCHEMES 



5 



The are nonnegative integers. On the other hand, since each E% o Ej (being a 
principal submatrix of Ei ® Ej) is positive semidcfinite, it follows that the are 
real and nonnegative. This important restriction is known as the Krein condition. 
The change-of-basis matrices P and Q are defined by 

n . n 

3=0 1 1 j=0 

We shall refer to P and Q as the first and second eigenmatrix of (X, 7?.), respectively. 
Note that Poi, Pu, . . . , P n i give the eigenvalues of Aj. The matrix P is also called 
the character table of (X, 7?.) , and in fact it can be viewed as a natural generalization 
of the character table of a finite groupH see Example 13.31 Let 

ki = P i, mi = Qoi- (4) 

It follows that hi is the valency of the regular graph (X, Ri) and ra; = trace(-Ej) = 
rank(£i). The rm are called the multiplicities of (X, 72). For convenience, set 
Afe = diag(fco, fci, • • • j fcn) and A m = diag(?no, mi, . . . , m n ). Then we have 

A m P = Q T A fe . (5) 

This is verified by evaluating tra,ce(AiEj) in two ways. Combining this with the 
obvious equality PQ = QP = \X\I, we get the following orthogonality relations: 

P r A m P=\X\A k , Q T A k Q = \X\A m . (6) 

We record here the eigenmatrix P for Examples O and O See [52 ] 155 1 155 1 120 ] [T4T] 
for the details. 

Example 2.3. Suppose (X,72) = J(v,n). Then fc, = ("X^"")' ^ = (*) - (.\) 
and the are given by dual Hahn polynomials [HFJ §1-6]: 

p ij _ P f-i,-j,i-v- 1 



fcj V n - u, — n 



Example 2.4. Suppose (X,7^) = H(n,q). Then fc, = m l = (")(<? - 1) J and the P- 
are given by Krawtchouk polynomials [951 §1.10]: 



p ij j? ( -i, -J 



—n 



q 



Wc remark that the polynomials in Examples 12.31 and 12.41 belong to the so-called 
Askey-schcmc of (basic) hypcrgcomctric orthogonal polynomials [95] , and the equa- 
tions in ([6]) for the corresponding association schemes amount to the orthogonality 
relations of these polynomials and their duals; see 33] 

An association scheme {X, S) on the same vertex set X is called a fusion of 
(X, TV) if each S € S is the union of some of the As the adjacency matrices 
(resp. primitive idempotents) of a fusion of (X, 1Z) must be 01-linear combinations 
of the Ai (resp. Ei), it is theoretically possible to find all fusions of (X, 72.) from 



°It is a long-standing problem whether the Pij are contained in a cyclotomic number field or 
not in general 1201 p. 123]. It is known that this is the case if the are rational ) 1301 ■ A negative 
answer to the problem would imply that the character theory of commutative association schemes 
is "far" from that of finite groups. See also 1961 . 

^By (AS1), -Ro £ S. As an extreme case, we call {X, IV) amorphous (or amorphic) if every 
"merging" operation on {Ri, R2, ■ ■ ■ , Rn} yields a fusion; see 1501 for a recent survey on this topic. 
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the eigenmatrix P. This is accomplished using the Bannai-Muzychuk Criterion 
[12llT34]H 



We close the section with a brief review of subschemes and quotient schemes. 
For a subset Y C X, define TZ Y = {Rj : < i «S n, Rj 7^ 0} where we write 
R Y = R n (Y x F) for R C X x X. We call (F,ft y ) a subscheme of (X,7l) if 
it is an association scheme. For example, J{v, n) can be naturally viewed as a 
subscheme of H(v, 2). 

We say that (X, 1Z) is primitive if the graphs (X, (1 ^ i ^ n) are connected, 
and imprimitive otherwise. Let I r (resp. J r ) denote the r x r identity (resp. all 
ones) matrix. Then 

Lemma 2.5. The following are equivalent: 

(i) (X, 7Z) is imprimitive. 

(ii) There is a subset T C {0, 1, . . . , n} such that X^ez Ai = I r ® J s for some 
integers r, s Jj 2 and an ordering of xV^ 

(iii) There is a subset J C {0, 1, . . . , n} smc/i i/ia£ X^ej ^« = s I r <8 J s / or 
some integers r, s Js 2 and an ordering of X. 

(iv) There exist i G {1, 2, . . . , n} and x,y € X (x ^ y) such that Eix — Eiy. 

(To prove this, proceed e.g., in the order (i)<^(ii)o(iii)=>(iv)=^(ii).) Suppose 
now that (X, TV) is imprimitive and let a subset X be as in (ii) above, so that 
lj ig 2: i?j is an equivalence relation on X . Each equivalence class affords a subscheme 
of (X, TV). We also have a natural structure of an association scheme on the set of 
all equivalence classes [501 §2-9], called a quotient (or factor) scheme of (X, 7£)q 
In fact, there is a concept of a "composition factor" of (X, 1Z) as well as a "Jordan- 
Holder" theorem. See [1771 1178] for the structure theory of (general) association 
schemes. Finally, it should be mentioned that if (X, 1Z) is primitive and symmetric 
then X can be viewed as a set of points on a sphere in M. mi for each i £ {1,2, . . . ,n} 
in view of (iv) above. This "spherical embedding" is known to be quite useful; see 
e.g., [30J Chapter 3], pH]. 

3. Gelfand PAIRS 

Suppose that a finite group G acts transitively on a finite set X (from the left). 
Let 1Z = {Rq, Ri, ■ ■ ■ , R n } be the set of the orbits^ of G on X x X under compo- 
nentwise action, where Rq = {(x,x) : x G X}, then (X,TZ) forms an association 
scheme Let n : G — > C XxX be the permutation representation of G on X, i.e., 
^{g)xy = S x ,g V (x,y & X, j e G) , and observe that the Bose-Mesner algebra A of 
(X, 1Z) coincides with the centralizer (or Hecke) algebra of 7r: 

A = {M G C XxX : n(g)M = Mir(g) for all g G G}. (7) 

Hence, by Schur's lemma it follows that (X, 1Z) is commutative if and only if 7r is 
multiplicity- free, i.e., if and only if tt is equivalent to a direct sum of incquivalent 
irreducible representations of G. We note that (X, 1Z) is symmetric if and only if 



^Scc e.g., 12 156 66, 83 for examples of explicit constructions of fusions using this criterion. 
^Such a subset is often said to be closed; see 11771 . 

12 Let E = s" 1 /,, @D J s be as in (iii) above. Then the Bose— Mesner algebra of the quotient 
scheme is canonically isomorphic to the "Hecke algebra" EAE (which is also an ideal of A). 
"^These orbits are also referred to as the orbitals or the 2- orbits of G on X ; see 1641 11751 . 
^Association schemes of this type are sometimes called Schurian 1641 . 
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7T is multiplicity-free and each irreducible constituent is realizable in R; this is also 
equivalent to the condition that the action of G on X is generously transitive, i.e., 
for any distinct x, y € X there is an element g £ G such that gx = y and gy = x. 

The G-set X can be identified with the set of left cosets G/K for some subgroup 
K of G. Note that n = (1k) g , where Ik is the trivial representation of K. For 
brevity, in this case we shall also let the symbol G/K denote the corresponding 
association scheme. It follows that G/K is primitive if and only if K is maximal 
in G. The pair (G, K) is called a Gelfand pair if (Ik) is multiplicity-free [1071 
Chapter VII], i.e., if G/K is commutative. If (G,K) is a Gelfand pair, then the 
determination of the eigenmatrices of G/K is in fact equivalent to describing the 
zonal spherical functions of (G,K) [20"1, §2.11]. 

The Johnson and Hamming schemes are obtained by actions of symmetric groups 
and their wreath products, respectively: 

Example 3.1. J(v,n) = S v /(S n x S v - n ). 

Example 3.2. H(n,q) = (S q I S n )/(S g -i I S n ). 

Concerning finite groups, we introduce another important example: 

Example 3.3 (The group association schemcF^I) . Suppose that X is a finite group 
and let G = Xxlactonlby (x,y)z = xzy~ l . Thus (X,H) = (XxX)/X in the 
above sense, where X — {(x,x) : x £ X}. We have (x,y) £ Ri if yx~ x € Gj, where 
Go = {1}, Gi, . . . , G n are the conjugacy classes of X. In this case A is isomorphic 
to the center of C x (viewed naturally as the group algebra of X) by At i— > J2 x ec 
Hence (X x X)/X is commutative and the Ei are in bijection with the irreducible 
characters £, ; of X. Let T be the group character table of X in the usual sense, and 
set A fc = diag(fc , ki, . . . , k n ), A/ = diag(/b, fi, ■ ■ ■ , f n ), where h = \d\ and f { is 
the degree of £j (0 ^ i ^ n). Then we havqj 

A f P = TA k , Q = T J A f . 

Note that (X x X)/X is primitive if and only if X is a simple group. 

The eigenmatrices of commutative association schemes have been extensively 
studied in the context of spherical functions. They are of significant interest in 
the theory of orthogonal polynomials (see e.g., |146j ). See also [1621,131?] for a wide 
variety of applications^ Besides, there are many other commutative association 
schemes whose eigenmatrices possess quite beautiful structure0 As a typical ex- 
ample we consider the association schemes 02 m +i(l) / 0^ m (q) , where for brevity we 



°E. Bannai (personal communication) pointed out that the concept of supercharacters and 
superclasses of finite groups introduced in 1601 is equivalent to that of fusions of group association 
schemes (in view of the Bannai-Muzychuk Criterion). 

^The central primitive idempotent corresponding to is given by fi\X\~^ JZoigx £ i{ x )^ see 
e.g., [48]. 

^The goal of the book 11621 is to consider finite analogues of the symmetric spaces including R" 
and the Poincare upper half plane, partly in order to "develop an understanding of the continuous 
theory by developing the finite model" 11621 p. 1], 

18 Concerning the character theory of finite groups, we especially refer to 12 1 1 1771 for a 
strong analogy between the character tables of GL(n,q), GU(n, q 2 ) and the eigenmatrices of 
GL(2n,q)/Sp(2n,q), GL(n,q 2 )/GL(n,q), GL{n,q 2 )/GU(n,q 2 ). See also [11]. 
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assume that q is evenF^ We first recall that 03(g) /0 2 (q) is symmetric with q/2 — 1 
classes (cf. [155] 179 ] ). The first eigenmatrix P is of the following form W\ 

9 + 1] 



where Pq is a square matrix of size q/2 — 1. (Recall that the top row of P 
gives the valencies.) Next, let e e {+,— } an d suppose m > 1 if e = — . Then 
02m+i(q) / 02 m (q) is symmetric with q/2 classes and the nontrivial valencies are 
given by h = (q^ 1 + (el))(q m - (el)), k 2 = k 3 = ■ ■ ■ = k q/2 = q m ~\q m - (el)). 
Moreover it turns out that the first eigenmatrix of 02m+i(q) /02 m (q) contains 
— (eg m_1 )Po as its lower-right submatrixri In this situation, we say that the eigen- 
matrix of 02m+i(<z) / 'O^miQ) * s controlled by that of the smaller association scheme 
Oz(q) / O2 (q) . There are a lot of fascinating examples of this kind; see [TT jll45| and 
the references therein. 

We close this section with an application of this theory. A connected /c-regular 
graph is called Ramanujan if all eigenvalues 9 such that \6\ 7^ k satisfy \9\ ^ 
2^/k — 1. Ramanujan graphs are good expanders and there are broad applications 
in computer science. Moreover, these graphs are precisely the regular graphs whose 
Ihara zeta functions satisfy the Ricmann hypothesis. See |162j for the details. 

Observe that 0^(q)/02~(q) = PGL(2, q)/D q+ i where D q+1 is a dihedral sub- 
group of order 2(q + 1), so that it is a quotient scheme of GL(2, q)/GL(l,q 2 ). The 
latter association scheme (for both even and odd q) is symmetric with q — 1 classes, 
and its q — 2 graphs with valency q + 1 arc called the finite upper half plane graphs 
[162] . These graphs are shown to be RamanuianF^ from which it follows that the 
(q + l)-rcgular graphs attached to O3 (g)/0^~ (q) are Ramanujan as well. An impli- 
cation of the above comments on the eigenmatriccs is that all graphs with valency 
q m ^ 1 (q m + 1) associated with 02m+i(9)/O^ m (<7) are Ramanuianl 23 ! see [23] . 

This systematic construction of Ramanujan graphs is an interesting application 
of the approach from the eigenmatrices of commutative association schemes^ In 
fact, the same method works for many other examples of controlling association 
schemes; see [22j [23]. See also [102[ l5T] for related constructions of Ramanujan 
graphs and [174j for an application of the results in [22l [23] to the Erdos distance 
problem. 



1 9 + 1 
1 



P 



Po 



^These association schemes arise from the action of 02m+l(g) on each of the sets of plus-type 
and minus-type hyperplanes. See |19l for the preceding results in the case of odd q. 

2 0ft is also known that all nontrivial multiplicities coincide (and equal q + 1), i.e., 03(g) / (q) 
is pseudocyclic 1301 §2.2B]. Pseudocyclic association schemes can be used to construct strongly 
regular graphs as well as distance-regular graphs with diameter three; see I3UI §12.7]. 

21 The proof is based on comparing the intersection numbers of these association schemes. Note 
that the other entries are determined from the orthogonality relations JB). 

22 The proof amounts to estimating several exponential sums over finite fields, such as Soto- 
Andrade sums, based on the work of A. Weil, N. Katz, W. Li and many others. See 11621 . 

2 ^Using a simple number-theoretic argument, it is also shown that for each fixed q there are 
infinitely many values of m such that the graphs attached to 02m+i ('?)/0^ m (g) are Ramanujan. 

24 Note, however, that this construction yields only finitely many Ramanujan graphs for each 
fixed valency. 
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4. METRIC/COMETRIC ASSOCIATION SCHEMES 

Suppose that (X, TV) is a symmetric association scheme. We say that (X, 1Z) is 
metric (or P- polynomial) with respect to the ordering {^4i}"=o ^ f° r eacrl i (0 ^ 
i ^ n) there is a polynomial Wj with degree i such that Pjj = Vi{Pj\) (0 ^ j ' ^ n). 
Such an ordering is called a P-polynomial ordering. Dually, we say that (X, 1Z) 
is cometric (or Q-polynomial) with respect to the ordering {Pi}" =0 if for each i 
(0 ^ i ^ n) there is a polynomial v* with degree i such that Qj; = v*(Qji) 
(0 ^ j ^ n). Such an ordering is called a Q-polynomial ordering. Note that in each 
of the above definitions the Vi (resp. v*) form a system of orthogonal polynomials 
by ([6]). Note also that (X,7V) is metric (resp. cometric) with respect to the above 
ordering if and only if for all i, j, k (0 ^ i, j, k ^ n) we have pfj = (resp. = 0) 
if i + j > k and ^ (resp. ^ 0) if i + j = k. 

A connected undirected graph (X, R) with diameter n and path-length distance 
d is called distance-regular if the n+ 1 relations Ri = {(.t, y) X x X : d(x, y) = i} 
(0 ^ i ^ n) define an association scheme [20l|30]. Thus metric association schemes, 
with specified P-polynomial ordering, are in bijection with distance-regular graphs. 
We refer the reader to [3D] for the basic theory on this topic0 and briefly comment 
on metric association schemes which are also cometric. (Henceforth we will use 
the phrase "metric & cometric" to describe such association schemes.) This class 
contains J(i),n), H(n,q) and many other important examples which arise from 
finite classical groups and classical forms over finite fields, e.g., Grassmann schemes 
and bilinear forms schcmes0 see J20J §3.6], (30J Chapter 9] and [49]. The famous 
theorem of Leonard [100] states that in this case the above polynomials Vi and v* 
belong to the terminating branch of the Askey scheme |95j (up to normalization), 
so that they are g-Racah polynomials [95j §3.2] or their special/limiting cases!3 see 
also [20l §3.5]. Recently, Leonard's theorem has been reformulated in the purely 
linear algebraic framework of Leonard pairs [1661 Q69]. Leonard pairs are used to 
describe certain irreducible modules for the Terwilliger algebra (Q of metric & 
cometric association schemes@ We also remark that Leonard pairs arise naturally 
in other various contexts, such as in representation theory^ see [Ml [HE] for the 
details. 

Compared with metric association schemes, the systematic study of cometric 
(but not necessarily metric) association schemes has begun rather recently. As we 
shall discuss below, they are of particular interest because of their connections, 
e.g., to spherical designs, (Euclidean) lattices and also mutually unbiased bases in 
quantum information theory. See also [17] . 



2 ^See also "Additions and corrections" to the book I30| available at A. E. Brouwer's webpage: 
|http:/ /www.win.tue.nl/~aeb/| (cf. footnote [3} . 

^These are g-analogues of J(v,n) and H(n,q), respectively. 

27 Note that we also allow the specialization q — » —1. Leonard's theorem was obtained just 
shortly after the 1979 discovery of the g-Racah (or Askey- Wilson) polynomials, and led Andrews 
and Askey [3] to their definition of the classical orthogonal polynomials. Bannai 1141 p. 27] 
states that "it is very interesting that combinatorics gave a meaningful influence to the theory of 
orthogonal polynomials in this way." 

28 Some problems on metric & cometric association schemes can be unified and most elegantly 
stated in terms of Leonard pairs; sec 159 . 



29 For example, we may obtain Leonard pairs from the finite dimensional irreducible modules 
for the Lie algebra sta as well as the quantum algebra U q (sl2)- 
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The famous Bannai~Ito Conjecture [201 p. 237] states that there are only finitely 
many distance-regular graphs with any given valency k > 2 (the polygons are all 
distance- regular with k = 2V For recent activity on this conjecture, see e.g., [571 [TU] 
and the references thercin|fj The dual to this conjecture is the following: 

Theorem 4.1 ( |123j ). For each fixed m > 2, there are only finitely many cometric 
association schemes with Q -polynomial ordering {£"i}" =0 satisfying rank(i?i) = in. 

A key step in the proof of this theorem is to bound the degree of the splitting 
field, based on the results of [150] ; see also [301 §3]. The splitting field of (X,TZ) is 
the smallest extension of the rational number field Q which contains all the Pij 

E 

While most distance-regular graphs with classical parameters have rational splitting 
field, the regular n-gon has splitting field Q(C) where £ = cos(27r/n) which, when 
n is prime for example, gives a degree (n — l)/2 extension of Q. In the case of 
distance-regular graphs of valency k > 2, only one known example — the Biggs- 
Smith graph has an eigenvalue not belonging to a quadratic extension of Q. To 
push this a bit further, to our knowledge, the only distance-regular graphs known 
with splitting field not contained in a quadratic extension of Q are 

• the Biggs-Smith graph, with spectrum 

31, 2 1 *, 17 , (±±^) 9 , (-l + 2cosf) 16 , 

(-1 - cos f + V3sin f ) 16 , (-f - cos f - V3 sin f ) 16 

• incidence graphs of generalized hexagons GH(q,q) 7 with eigenvalues 

±(g+l), 0, ±^q, ±y/3q. 
(These can be viewed as point graphs of thin generalized 12-gons of order 

• line graphs of the above graphs, with eigenvalues 

2q, q-1, -2, q-l±y/q, q - 1 ± yfa, 
which are the point graphs of generalized 12-gons of order (g, 1). 
The question arises as to whether there exists a bound on the dimension of the 
splitting field of a distance-regular graph in terms of its valency. Any such result 
would immediately imply the Bannai-Ito Conjecture since the interval [—k, k] would 
then contain a limited number of potential eigenvalues (since it must contain all 
their conjugates). But such a bound seems beyond our reach at this point. In [30[ 
p. 130], Brouwer, et al. ask if a distance-regular graph with k > 2 must have at least 
two integral eigenvalues. Even this apparently simple question remains unresolved 
to date. 

Conjecture 4.2 (Bannai and Ito [201 p. 312]). For n sufficiently large, a primitive 
association scheme with n classes is metric if and only if it is cometric. 

While no counterexamples are known to this conjecture, there are imprimitive 
distance-regular graphs — the doubled Odd graphs — which are not cometric yet 
have arbitrarily large diameter, and there are some Q-bipartite doubles of certain 
Hermitian forms dual polar spaces [20i p. 315] which are cometric yet not metric. 



30 As of this writing, it has been announced that the full conjecture has been proven by Bang, 
Koolen and Moulton. 
31 See footnote [8] 
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So the "primitive" condition in the conjecture is necessary. Still, no progress has 
been made toward proving the conjecture. From the "cometric viewpoint," perhaps 
the following questions will prove easier to attack: 

Question: Suppose that (X, 1Z) is cometric with Q-polynomial ordering {Ei}f_ 
and the A, ordered so that Qoi > Qu > • • • > Qni- If n is sufficiently large, must 
the adjacency matrix A\ of relation R\ have n + 1 distinct eigenvalues? 

Question: Does there exist an absolute constant v such that, for any cometric 
association scheme (X, TV) with Q-polynomial ordering {Ei}™ =0 and the Ai ordered 
so that Qoi > Qn > ■ ■ ■ > Q n i, we have p\- = whenever k > j + vl 

It is easy to prove that the valencies of any distance-regular graph with diameter 
n satisfy the unimodal property: 

k = 1 < fei 5$ k 2 ^ • • • s% fe„' ^ fen'+l > . • • > fen (8) 

for some 1 ^ n' ^ n, possibly n' = n. For cometric association schemes, we have 

Conjecture 4.3 (Bannai and Ito [20l p. 205]). If (X,TZ) is cometric with Q- 
polynomial ordering {Ei}™ =0 , then the mi satisfy the unimodal property: 

mo = 1 < mi ^ 7712 ^ • ■ ■ ^ mn' 3 s fl^n'+l 3 s ■ • ■ 3 s m n 

for some 1 ^ n' ^ n. 

In personal communication with P. Terwilliger, the following stronger claim was 
made for all cometric association schemes: 

Conjecture 4.4 (D. Stanton). For i < n/2, mi ^ and m.; ^ m n -i. 

Stanton's conjecture has been proven under the added assumption that the as- 
sociation scheme is dual thin (^J9| |140j or metric [138] . 

It is well known that a metric association scheme admits at most two P-polyn- 
mial ordcrings [30, Theorem 4.2.12]. §4. 2D in [301 examines the possbilitics for a 
second P-polynomial ordering of a distance-regular graph and obtains substantial 
parameter conditions. In the cometric case, we have the following result of Suzuki 
(1998): 

Theorem 4.5 ( |150| ). If {Ei}f_ is a Q-polynomial ordering for a cometric asso- 
ciation scheme (X,TZ), then any second such ordering must be one of: 

• E , E 2 , E4, . . . , E 3 , Ei; 

• Eq, E n , Ei, E n _i, . . .; 

• Eo,E n ,E 2 ,E n - 2 ,--- (n odd); 

• Eq, E n -i, E 2 , E n -3, ... (n even); 

• Eq,E 5 ,E 3 ,E 2 ,E4,Ei. 

Further conditions were given in |150j (e.g., (X, 1Z) must be almost Q-bipartiteF^ 
in the first case). It is possible that the last case may be ruled out. 

The Krein parameters of a cometric scheme (X, 7Z) are entirely determined by 
its Krein array 

l*(X,TZ) = {b* ,b*i, . . . ,b n _i;cl,c 2 , . . . ,c* n }, (9) 



A cometric association scheme with n classes is almost Q -bipartite if its parameters satisfy 
a* = for all i < n, yet a* > 0. 
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where b* = q\ i+ i (0 ^ i ^ n — 1) and c* = q\ i _ l (1 ^ i ^ ri). We also define 
a* = q l u (0 ^5 i $J n). It is well known [20) p. 315] that an imprimitivc distance- 
regular graph with valency k > 2 is bipartite or antipodal (or both). The dual 
situation is not yet fully resolved. In 1998, Suzuki proved 

Theorem 4.6 f [149j ). If (X,1Z) is an imprimitive cometric association scheme 
with Q -polynomial ordering {Ei}" =0 , then at least one of the following holdsW^ 

(i) (X, 1Z) is Q-bipartite: a* = for 1 ^ i ^ n; 

(ii) (X,1Z) is Q-antipodal: b* ~ c* —i for 1 ^ i ^ n, except possibly i = [n/2\; 

(iii) n = 4 and l* (X, 1Z) = {to, to — 1, 1, 63; 1, Cj, to — 63, 1}, where a 2 > 0; 

(iv) n = 6 and i*(X, TV) = {m, m — 1, 1, 63, b\, 1; 1, c^, m — 63, 1, C5, to}, where 
al = a'l + al > 0. 

Schemes of type (iii) in the theorem have recently been ruled out [3D]. No 
examples are known of type (iv); it is quite possible that none exist and then the 
theorem exactly mirrors the result for imprimitive distance-regular graphs. 

Let us briefly review the known examples of such imprimitivc "polynomial schemes' 
with three or four classes. In the metric case, a bipartite distance-regular graph of 
diameter three is necessarily the incidence graph of some symmetric (v, k, A) block 
design. Any such scheme is cometric as well. An antipodal distance-regular graph 
of diameter three is a cover of a complete graph [74) ; these are cometric precisely 
when the cover has index two. The only distance-regular graphs of diameter three 
which are both bipartite and antipodal are the complete bipartite graphs with a 
perfect matching deleted. These trivial examples arc both metric and cometric. 

The 3-class imprimitivc cometric schemes follow a landscape dual to this. The 
Q-bipartite examples are all Taylor graphs; they are all index two distance-regular 
antipodal covers of the complete graphs. A 3-class Q-antipodal scheme is equivalent 
to a linked system of symmetric designs |361 1119) ; these are only metric when there 
are two Q-antipodal classes, these being the incidence graphs of symmetric designs 
mentioned above. The only examples which are both Q-bipartite and Q-antipodal 
are again the complete bipartite graphs with a one-factor deleted. 

In the case of imprimitive 4-class schemes, the bipartite distance-regular graphs 
of diameter four are incidence graphs of various designs and geometries (e.g., gen- 
eralized quadrangles) while the antipodal distance-regular graphs of diameter four 
are antipodal covers of strongly regular graphs. The distance-regular graphs of di- 
ameter four which are both bipartite and antipodal arc characterized as incidence 
graphs of symmetric (to, /i)-nets [301 P- 18] ■ On the cometric side the 4-class schemes 
which arc Q-bipartite correspond to systems of lines with two angles, one of which 
is 7r/2; the 4-class Q-antipodal schemes are roughly the linked systems of strongly 
regular designs. Interestingly, the 4-class schemes which are both Q-bipartite and 
Q-antipodal are in one-to-one correspondence with sets of real mutually unbiased 
bases [HFJ, which we now define. (See also [T].) 

Let {61, 62, ... , bd} and {b[, b' 2l . . . , b' d } be two orthonormal bases for C d . We 
say these bases are unbiased (relative to one another) if |(6j|6j-)| = l/Vd for all 
1 ^ i, j ^ d where (|) is the standard Hermitian inner product on C d . A collection 
of orthonormal bases for C d is mutually unbiased if any two distinct bases from the 
set are unbiased relative to one another. For d a prime power, there is a construction 



■^With the notation of Lemma 12.51 the types (i)-(iv) correspond to J = {0, 2, 4, . . . }, J = 
{0, n}, J = {0, 3} and J = {0, 3, 6}, respectively. 
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of d + 1 mutually unbiased bases (MUBs) in C d . For other dimensions this is 
mostly an open question. Such constructions are useful for several applications in 
quantum information theory, such as quantum key distribution and quantum state 
tomography. (See [29] and the references therein.) 

When we restrict the bases to be real, the problem changes qualitatively. Indeed, 
for unit vectors b and 6' from distinct bases, we must have {b\b') — ±l/yd. Let 
Md denote the maximum possible number of real MUBs in dimension d. It follows 
from an old result of Delsarte et al. [55] that Md ^ d/2 + 1; see also [32]. This 
bound is achieved for d = 4 fc via a construction using Kerdock sets. In [25], it is 
established for example that 

• Md = 1 unless d = 2 or 4\d; 

• Md ^ 2 if an only if there exists a Hadamard matrix of side d; 

• Md < 3 unless d/4 is an even square. 

Using the results of [99] , each of these results gives either a construction or a non- 
existence result for cometric association schemes with four classes which are both 
Q-bipartite and Q-antipodal; specifically, Md is an upper bound on the number k 
of Q-antipodal classes in such association scheme on 2kd vertices with Q-antipodal 
classes of size 2d. 

If (X, 1Z) is Q-bipartite with Q-polynomial ordering {Ei}2_ Q , then the set {E\x : 
x G X} is closed under multiplication by — 1; so, viewed as points on the unit 
sphere in M. mi , these schemes are really best viewed as sets of lines through the 
origin. The imprimitivity system here has all equivalence classes of size two and 
the quotient scheme, on |X|/2 vertices, is often interesting. Examples include the 
schemes arising from the shortest vectors in the Eq, Et, E$ and Leech lattices, as 
well as an overlattice of the Barnes- Wall lattice in IR 16 ; these have Krein arrays 

. l*(E 6 ) = {6,5,1,3; 1,§,3,6} 

. ,*(£ 7 ) = { 7 ,6,f,f;l,f,f,7} 

. l*(E 8 ) = {8,7, f ,6; 1,|,2,8} 

. S (Leech) = {24, 23, m IfO, 104, si. ^ * ^ M ^ 24 } 

• t*(OBW16) = {16,15, 8; 1,^,8, 16} 

Further Q-bipartite examples come from Q-bipartite doubles of certain strongly 
regular graphs, such as the two subconstituents of the McLaughlin graph. 

Concerning the structure of Q-antipodal schemes, again very little is known. 
The quotient scheme is a one-class scheme. With the natural ordering on the Ai, 
we have T = {0, 2, 4, . . . } in the notation of Lemma |2.5l The following theorem has 
been referred to as the "Dismantlability Theorem" ri 

Theorem 4.7 f |119| ). If (X,7Z) is Q-antipodal and Y C X is a union of £ Q- 
antipodal classes, then (Y,TZ Y ) is a cometric subscheme, which is Q-antipodal as 
well, provided £ > 1 . 

In fJB] we shall investigate duality among association schemes. As a special 
case, if (X,TZ) is the coset scheme of an additive completely regular (§5} code Y 
in H(n,q), then the dual of (X,TZ) (induced on the dual code Y°) is a cometric 
subscheme inside H(n,q). In this way, we obtain a number of cometric schemes 



^ This theorem is formally dual to an unpublished result of C. Godsil (personal communication) 
which states that, in an antipodal distance-regular graph, any subset of an antipodal class is a 
completely regular code (^S). 
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from the perfect binary and ternary Golay codes and some codes derived from 
them [30j p. 356]. What is new here is that, since several of these coset graphs are 
antipodal, their dual schemes are Q-antipodal anc j the above theorem gives us new 
cometric schemes which are not metric. We give two examples here. 

• The dual scheme of the coset graph of the shortened ternary Golay code is 
a Q-antipodal scheme on 243 vertices with three Q-antipodal classes. If we 
dismantle this, taking two of these classes only, we obtain a Q" an t.ipodal 
scheme with Krein array {20, 18, 3, 1; 1, 3, 18, 20}. It is interesting to note 
that the dual parameter set remains open for a possible antipodal diameter 
four distance-regular graph. 

• Example (A16) on p. 365 of [30] is the coset graph of an additive binary 
code derived from the extended binary Golay code. Its dual scheme has 
2048 vertices and four Q-antipodal classes. If we take only three of these, 
we obtain a Q-antipodal scheme with Krein array {21, 20, 16, 8, 2, 1; 1, 2, 4, 
16,20,21}. In this case, the dual parameter set has been shown to be 
unrealizable as a distance-regular graph by counting hexagons in such a 
graph .'!() p. 365]. 

In a terse summary of spherical designs [13 1J . Munemasa gives numerous ex- 
amples of cometric schemes arising from lattices which are not distance-regular 
graphs. Martin et al. [119] build on this list, including some schemes coming from 
error-correcting codes, block designs and the above theorem applied to known Q- 
antipodal schemes. Higman's paper [78] on strongly regular designs contains further 
examples. 

We have already mentioned some imprimitivc examples. It is remarkable that 
very few primitive cometric association schemes are known which are not metric. 
The only known examples, to our knowledge, are the following Fl* 

• the block scheme of the 4-(ll, 5, 1) Witt design, with n = 3, \X\ = 66 and 
Krein array {10,242/27,11/5; 1,55/27,44/5} 

• the block scheme of the 5-(24,8, 1) Witt design, with n = 3, \X\ = 729 
and Krein array {23, 945/44, 1587/80; 1, 345/176, 207/20} 

• a spherical design derived from the Leech lattice with n = 3, \X\ = 2025 
and Krein array {22,21,625/33; 1,11/6,30/11} 

• the block scheme of a 4-(47, 11,48) design arising from codewords of weight 
11 in a certain quadratic residue code of length 47, with n — 3, \X\ = 4324 
and Krein array {46, 77315/1782, 24863/847; 1, 37835/19602, 2162/231} 

• the "antipodal" quotient of the association scheme on shortest vectors of 
the Leech lattice, with n = 3, \X\ = 98280 and Krein array 

{299, 1800/7, 4563/20; 1, 156/35, 195/4} 

• a spherical design derived from the Leech lattice with n = 4, \X\ = 7128 
and Krein array {22,21,121/6,2187/125; 1,11/6,363/125,6} 

• another derived spherical design arising among the shortest vectors of the 
Leech lattice, with n = 5, \X\ = 47104 and Krein array 

{23, 22, 529/25, 184/9, 483/25; 1, 46/25, 23/9, 92/25, 23/3} 



"^See also an on-line table of cometric association schemes which are not metric, maintained 
by W. J. Martin: |http://users.wpi.edu/~martin/RESEARCH/QPOL/| 
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5. Codes and designs 

Suppose that (A, 1Z) is a commutative association scheme. Throughout this 
section, let Y be a nonempty subset of X with 1 < \Y\ < \X\. Let \ = Y^xeY^ 
be the characteristic vector of Y . The inner distribution of Y is the vector a = 
(ao, Oi, . . • , a n ) defined by 

a ' = ^ Ti ' x = M |fl ' n(yxy)l ' (10) 

Note that the a^ are nonnegative, ao = 1 and {aQ)o = \Y\^ Since (aQ)i = 
|A| \Y\~ 1 x T EiX it follows that the (aQ)i are also real and nonnegative; this simple 
fact underlies Delsarte's linear programming method; see SJ7] We remark that 
(aQ)t = if and only if Ei\ = 0. 

For a subset T of {1, 2, . . . , n}, we call Y a T-code (resp. (Delsarte) T-design) 
if a; = (resp. (aQ)i = 0) for all i G T. A {1, 2, ... , i}-design is simply called a 
t-design\^\ Codes in H(n, q) are the familiar "block codes of length n," and codes 
in J(v, n) are precisely the binary constant-weight codes. We remark that codes in 
the bilinear forms schemes also have applications to space-time codes; see |76j . 

If Y is a T-code and if Z C X is a U-code with inner distribution b where 
T U U = {1, 2, . . . , n}, then by the right side of (O we have 

\Y\\Z\ < (aQ)A~ 1 (6(3) T = |A>A^6 T = LY| (11) 

with equality if and only if (aQ)i(bQ)i = (1 ^ i ^ n). This "Anticodc Bound" 
is a special case of the linear programming method. A similar argument gives an 
"Antidesign Bound" for T-designs: if Y is a T-dcsign and Z is a W-dcsign where 
TUW={l,2,...,n}, then |F||Z| ^ |A|. 

In some cases, certain T-dcsigns have natural geometric interpretations. For 
example, if (A, 7?.) is induced on the top fiber of a shoriF^ regular semilattice ("P, 
(see [13]), then Y is a i-desigrj3 if and only if the number € Y : u =^ x}\ (called 
the index) is independent of u G V with rank(it) = t [53] . For Examples 12.11 and 
12.21 we have: 

Example 5.1. Let V = {u C f2 : |u| ^ n}. Then (V,^), where the partial order 
is given by inclusion, forms a short regular semilattice (truncated Boolean lattice) 
with rank function rank(u) = \u\. In the top fiber J(v,n), a Delsarte f-design is 
just a combinatorial i-design 

Example 5.2. Introduce a new symbol "■" and let V be the set of words of length n 
over Q U {•}. For u = (u\, U2, ■ ■ ■ , w n ), u = (v\,V2, ■ ■ ■ , v n ) G T 3 , we set w, =<l i; if and 
only if Ui = ■ or Ui — z>i, for all i. Then (V, =4) defines a short regular semilattice 



In general, for a vector c = (co, ci, . . . , c n ) we call cQ the MacWilliams transform of c. 
^ In what follows, if we define a concept /parameter which depends on the ordering of the Ai 
or the Ei (such as a t-dcsign) then we shall understand that such an ordering is implicitly fixed. 
Whenever we state a result involving these concepts/parameters, the orderings will be explicitly 
specified or clear from the context. 

ranked, meet semilattice (V, =4) with top fiber X is short if X A X = V . 
■^Here we are using the ordering defined naturally by the semilattice structure. 

^"A t-(v, n,A) design is a collection of n-subsets (called blocks) of a v-sct such that every 
t-subsct is contained in exactly A blocks. 



16 



WILLIAM J. MARTIN AND HAJIME TANAKA 



(Hamming lattice) with rank function rank(u) = \{i : m ^ -}|. In the top fiber 
H(n, q), a i-design is an orthogonal array of strength t ¥^\ 

See |129[ I148| for geometric interpretations of i-designs in other classical families 
of metric & cometric association schemes. More "exotic" types of codes and designs 
will be discussed in Sj8] See also |33l [57] for another approach to the regularity of 
T-designs in J(v,n) and H(n,q) in terms of their t-form spaces. 

The outer distribution of Y is the \X\ x (n + 1) matrix 

B = [A oX ,A 1 x,...,A nX ]. (12) 

We also recall the following four fundamental parameters of Y : 

8 = mm{i ^ : a, ^ 0}, 8* = min{i ^ : (aQ) l ^ 0}, (13) 

s = \{i^0: ai ^0}\, 8* = \{i?0:(aQ)i?0}\. (14) 

We call 8, 8*,s,s* the minimum distance, dual distance, degree and dual degree of 
Y, respectively^ These are related with |Y| as follows: 

Theorem 5.3 f |52|). Suppose that (X,1Z) is metric with P -polynomial ordering 
{Ai}™ =0 . Then 8 < 2s* + 1 and 

L(5-l)/2j 

E fc ^U^E fc - ( 15 ) 

i=0 1 1 i=0 

If 6 ^ 2s* — 1 £/ien Y is completely regular, i.e., the X th row of B depends only on 
0(x,Y) = mm{i: B xt ^0}. 

Theorem 5.4 ([52]). Suppose that (X,1Z) is cometric with Q-polynomial ordering 
{Ei}f =0 . Then 8* ^ 2s + 1 and 

L(«5*-1)/2J 

53 m<<|y|<53m i . (16) 

i=0 i=0 

7/(5* ^ 2s — 1 i/iera (Y, 72. ) is a cometric subscheme with s classes. 

The inequality in the left side in (|15[) (resp. (|16p ) is the sphere-packing bound 
(resp. Fisher-type inequality), and Y is a perfect code (resp. design) if it 

satisfies equality. It follows that Y is a perfect code (resp. tight design) if and only 
if 8 = 2s* + 1 (resp. 8* = 2s + 1). We remark that the codes with 8 £ {2s* - 1, 2s*} 
in Theorem 15.31 are precisely the uniformly packed codes |30[ p. 348]. Completely 
regular codes have been actively studied because of their importance in the theory 
of distance-regular graphs; see Chapter 11] and [116] . 

Suppose now that (X,7Z) is metric with P-polynomial ordering {Ai}l l =Q . Pick 
any x 6 X and set Z x = {y £ X : (x,y) £ {Ji =0 Ri}, wncre e — 1(8 — 1)/2J is 
the packing radius of Y. We may obtain the sphere-packing bound via (jlip with 
Z = Z x . Thus, if Y is perfect, then since the characteristic vector ip x — "YT i=0 A^x 

of Z x satisfies ^ X E^ X = (YZ=oPii? m i\ x \~ 1 > we find \U ^ : E<=o^ = °>l = 
e(= s*). In other words, using the notation at the beginning of 21 all the zeros 
of the Lloyd polynomial y,?_ n i>j must be in {Pn, P21, . . . , P n i}- We remark that 
this "Lloyd Theorem" has a dual, so that we also obtain a strong nonexistence 

^An orthogonal array OA\{t,n,q) is a Ag' X n matrix over an alphabet Q of size q in which 
each set of t columns contains each t-tuple over Q exactly A times as a row. 
^ 2 We also refer to r = 8* — 1 as the (maximum) strength of Y. 
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condition on tight designs in general cometric association schemes in terms of the 
Wilson polynomial Yli=o v t- See for the details 

Wc may derive a lot more structural information on Y by just looking at the 
four parameters. For example, it follows from BQ — \X\[EoX, E1X1 ■ ■ ■ , E n x] that 
rank(P) = s* + 1. Hence, if (X,1Z) is metric with P-polynomial ordering {A;}" =0 , 
then the covering radius p = max{9(x, Y) : x G X} of Y must satisfy p ^ s*. (The 
right side of (fT5|) follows from this.) We call Y regular if x T AiX is independent of 
x G Y (and thus equals dj) for all i. It is known that 

Theorem 5.5 ([52]). Suppose that (X,1Z) is metric with P-polynomial ordering 
{Ai}™ =0 . If S s* then Y is regular. 

Theorem 5.6 ([52]). Suppose that (X,1Z) is cometric with Q-polynomial ordering 
{EiYl =0 . If 3* s then Y is regular. 

We refer the reader to |59[ I37] for more detailed information and the (major) 
progress up to 1998. We remark that Delsarte's theory of codes and designs (in 
metric/cometric association schemes) based on the linear programming method has 
been naturally extended to various compact symmetric spaces with rank one 



such 

as spheres S n = SO(n + l)/SO(n); see fTT] for a survey on this topic. 
In 2003, Brouwer et al. [31] introduced the following parameters for Y: 

w = max{i : a; ^ 0}, if* = max{i : {aQ)i ^ 0}. (17) 

We call w, w* the width and dual width of Y, respectively. They obtained the 
following results: 

Theorem 5.7 (|31j). Suppose that (X,1Z) is metric with P-polynomial ordering 
{AiYi = Q- Then w ^ n — s* . Ifw = n — s* then Y is completely regular. 

Theorem 5.8 ([31]). Suppose that [X,1Z) is cometric with Q-polynomial ordering 
{Ei}i=o- Then w* ^ n — s. If w* = n — s then (Y,TZ Y ) is a cometric subscheme 
with s classes. 

The above results are in contrast with the bounds on S and 5* in Theorems 
15.31 and 15.41 See [31j [80] for many interesting examples attaining the bounds in 
Theorems 15.71 and 15.81 

Suppose now that (X,1Z) is metric with P-polynomial ordering {j4i}™ =0 and 
cometric with Q-polynomial ordering {£"i}™ =0 . Since w s and w* ^ s* wc have 

w + w* ^ n. (18) 

If (X, TV) is induced on the top fiber of a short regular semilattice (V, =4), then for 
every u G V the subset Y u = {x G X : u =4 x} satisfies w = n — rank(u) and 
to* = rank(u)0 It is shown in [31] that any code with w + w* = n in J(y,n) 
and H (n, q) is isomorphic to a code of the form Y u . This result was later extended 
to their g-analogues in [157] , It should be remarked that for these examples the 



In fact, one may find an analogue of this Wilson polynomial in schemes which are not 
necessarily cometric. Let £, T C {0, 1, . . . , n}. Define £ * T to be the set of eigenspace indices k 
(0 ^ k $J n) such that g*. ^ for some i G £ and some j S T . Then, if Y is a T-dcsign and £ 
satisfies £ * £ C T U {0}, we obtain the lower bound \Y\ ^ 5Zjg£ m j- See 11151 for details and 
further conditions on the case when equality holds. 

^In this case, the corresponding orthogonal polynomials are Jacobi polynomials 1951 §1.8]. 

^The characteristic vectors of the Y u with rank(u) = £ span JZi=o Yi (0 ^ & ^ n); see |53j . 
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Y u again afford metric & cometric association schemes which belong to the same 
family as the orig inal0 This "hierarchical structure" appears to be a subject ripe 
for further investigation. 

The Erdos-Ko-Rado Theorem [63] states that for each integer t such that v > 
(t + l)(n — t + 1) the largest codes satisfying w ^ n — t in J(v, n) are the Y u with 
rank(w) = t. In fact, the original proof in [63] based on the "shifting technique" 
establishes the conclusion under the stronger hypothesis v ^ t+ (n — £)(") , and the 
best possible bound v > (t + l)(n — t+1) was obtained in [176j as an application 
of Delsarte's linear programming method. The observation that the largest (or 
extremal) codes in the Erdos-Ko-Rado Theorem are those codes satisfying w+w* = 



sly i 




n led to the "g-versions" of the theorem in full generality; see |157j 

6. Duality 

Suppose that (X, TV) is a commutative association scheme and that X is endowed 
with the structure of an abelian group (written multiplicatively) with identity el- 
ement 1. We call (X, TV) a translation association scheme if for all ^ i ^ n and 
z G X, {x,y) G Ri implies (xz,yz) G Ri. This concept is equivalent to that of a 
Schur ring on an abelian group; see [136j for a survey on Schur rings. 

Let X* be the character group of X. To each e G X* we associate the vector 
i = \X\- 1 / 2 J2 xe x^( x ~j x > so that (M) = \X\- 1/2 e(x)^ Note that the e form an 
orthonormal basis for C x by the orthogonality relations for the characters. Define a 
partition X = {X 0l X u ..., X n } of X by Xi = {x G X : (1, x) G Ri} (0 < i n)@ 
Then Ri = {(x, y) G X x X : yx- 1 G Xi} (0 i sC n) and we find 



Aii = [ } e{x) ]i (0 < i < n, e G X*). (19) 




Hence we may also partition X* as follows: X* = {Xq , X*, . . . , X*}, where X* = 
{s e X* : e G Vi} (0 sC i s$ n). It follows that 

= J2 ( £ e o« = E £ (- t ) (* e x ») ( 2 °) 

for ^ i,j ^ n. The left-hand equation of (|20[) is immediate from (HHJ), and 
the right-hand equation follows by evaluating |X|(-Ej)i.j; in two ways using Ej = 
J2eex* e^ T Let K* = R\ , . . . , R* n } be the partition of X* x X* defined by 

R* = {(e,V) : V^ 1 G (0 < i s$ n), and let A* G C A "* xX * be the adjacency 
matrix of (X*,R*) (0 sC i ^ n). If we identify C x ' xX ' with C XxX via the or- 
thonormal basis {e : e G X*}, then it follows from the orthogonality relations and 
AH that 

n 

A t= E ^ T = E^*' ( 21 ) 



4 ^At the algebraic level, this is explained from the results in 11591 . 

^This is a consequence of the previous work 1651 1811 , together with the classification of codes 
with w + w* = n. The construction of Singleton systems 1561 (i.e., t-designs with index one) in 
bilinear forms schemes plays an important role in the proof (in view of Him ): see also 1811 p. 192], 

48 See footnote [T6l 

^Such a partition is sometimes referred to as a blueprint; see e.g., [8]. 
^Note that {e : e £ X*} forms an orthonormal basis for Vi. 
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where E* is the diagonal matrix in <C XxX with (cc, x)-entry (E*) xx = (Ai)i x ; so 
the vector space A* spanned by the A* is an algebra. Hence (X* , TV ) is again a 
translation association scheme, called the dual of (X,TZ). By (f!H]) , (X*,1Z*) has 
cigcnmatriccs P* = Q and Q* = P. This duality was first formulated in |152j . but 
the structure of the Terwilliger algebra ($9| is already visible here. 

Let Y be a subgroup of X with characteristic vector x an d inner distribution 
a = (ao, eti, . . . , a„). Note that Y is regular and thus at = \Y D (0 i ^5 n). 
Set r = {eer: e(j/) = 1 for all y G Y}. Then 7° is a subgroup of X* and 

(aQ), = ^fx T ^X=|f| £ |<X,e*>| 2 = m-|l r °nJf;|. (22) 

v* 

It follows that Y° has inner distribution |Y| _1 a(3. We remark that H(n,q) is a 
translation association scheme if we take the alphabet Q to be an abelian group (cf. 
Example Moreover, in this case the dual of H(n, q) is again the Hamming 

scheme (with vertex set X* = (Q*) n ); in other words, H{n,q) is self-dual. Thus, 
in view of the generating functions for the Krawtchouk polynomials §1.10], 
(|2"2"| turns out to generalize the well-known MacWilliams identity on the weight 
distributions (or enumerators) of a linear code and its dual code. The following 
theorem is also important (cf. Theorems 15.41 and 15.81) : 

Theorem 6.1 (|52|). With the above notation, (Y,1Z Y ) is a subscheme if and only 
if the outer distribution ofY° has s + 1 distinct rows, where s is the degree ofY. 

If (Y, TZ Y ) is a subscheme (with s classes), then its dual scheme has vertex set 
X* I Y° and the relation containing a pair (eY°, rfY°) is determined by the (?7£ _1 ) th 
row of the outer distribution of Y° ; see [52] ■ 

Certain dualities of commutative (but not necessarily translation) association 
schemes also arise in connection with spin models and type II matrices. Let A be 
a nowhere zero matrix in C XxX with "Schur inverse" A^ , i.e., A o A^ = J. 
(Henceforth we shall not assume a group structure on X .) We call A type II if 
AA^ T — \X\I. The Nomura algebra of A is the space Ma of matrices M in C XxX 
such that Ax o A^y is an eigenvector of M for all x,y £ X. HA is invertible, 
then A is type II if and only if J G Ma (cf- [131 Lemma 2.1]). Define a linear map 
S A : M A -> C XxX by 

M(Ax o A^y) = (e A {M)) xy ■ {Ax o A^y) (M G M A , x, y G X). 

Jaeger ct al. [SD] showed that if A is type II then Qa{Ma) — M/p, and 

e A T{Q A (M)) = \X\M T 1 e A (MN) = e A (M)oQA(N) (23) 

for M, N G AOi ■ It follows that if A is type II then both Ma and A/^t are the 
Bose-Mesner algebras of some commutative association schemes, and ®a gives an 
isomorphism between them which "swaps" the ordinary multiplication and o. 

Spin models were introduced by V. Jones |91| as a tool for creating link invariants, 
and are characterized (up to scalar multiplication) as those type II matrices A 
satisfying A G Ma ESI Proposition 9]. If A is a spin model then in fact we have 
Ma = Mat and 6^t = 0^ [90l Theorem 11]; in this case ([23]) is equivalent to 
the condition that the corresponding association scheme is formally self-dual, i.e., 
P = Q for some orderings of the Ai and the see [18]. In fact, it was shown that 



The most familiar case is that Q is a finite field F 9 and Y is a linear code in the usual sense. 
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Ma is not just formally self-dual, but is "strongly hyper-self-dual" which is defined 
using the Terwilliger algebra; see [47]. See e.g., [38] [46] for more information on 
the connections to the Terwilliger algebra and Leonard pairs. Spin models, as well 
as four-weight spin models |15j . have been studied via the more general but crisp 
concept of Jones pairs [44] [44] [42] . 

A formally self-dual association scheme (X, TV) is said to satisfy the modular 
invariance property (with respect to P) if there is a diagonal matrix A such that 
(PA) 3 is a nonzero scalar matrix. This gives a necessary condition that A = Ma 
for a spin model A [181 190] . The modular invariance property is also quite relevant 
to fusion algebras in conformal field theory [13, 68]. We remark that fusion algebras 
are closely related to table algebras [25] and to character algebras [20] §2.5] which 
may in turn be viewed as "Bose-Mesner algebras at the algebraic level." 

7. The linear programming bound 

Suppose that (X, 1Z) is a symmetric association scheme. Let T be a subset of 
{1,2, ... ,n}. In coding theory, we are often interested in finding a sharp upper 
bound on the size of a T-code in X. The fact that the inner distribution a of 
a code and its "MacWilliams transform" aQ are nonnegative leads to the linear 
programming (or LP) bound developed by Delsarte: 

Theorem 7.1 ([52]). With variable a = (ao, oi, . . . , a n ) £ M. n+1 , set 

4p = 4p (X, T) = max (oQ) (24) 

subject to (i) oq = 1; (ii) a, ^ (1 ^5 i n); (Hi) (aQ)i (1 ^ i ^ n); (iv) 
ai = 0ifie T. IfY C X is a T-code, then \Y\ ^ £ L p- 

The LP bound was shown to be a close variant of Lovasz's $-bound |105| on 
the Shannon capacity of a graph; see |142j . Many computational software packages 
implement the simplex method to solve linear programming problems, and (12"4"]1 
does produce a lot of sharp upper bounds on the size of codes. However, most 
analytic results give bounds using the dual linear program: 

Theorem 7.2 ([52]). With variable b = (b 0) h, b n ) E W n+1 , set 

4p = 4p (X, T) = min (bQ T ) (25) 

subject to (i) b Q = I; (ii) b l > (1 sC i ^ n); (Hi) (5Q T ); < ifi G {1, 2, . . . , n}-T. 
Then £lp = ^lp ■ 

If a, b are feasible solutions to the programs ([2"4"]l and ([2"S"]) respectively, then 

(oQ) < aQb T ^ (bQ T ) (26) 

with equality if and only if (aQ)ibi = a.i(bQ T )i = for 1 ^ i ^ n. Note that 
([TT]) also amounts to constructing a feasible solution to ([25]) ; we saw in f]5] that 
the optimality condition was the key to prove and generalize Lloyd's Theorem. In 
passing, ([TT]) can be slightly strengthened as follows: i^{X, T)£i I p(X,U) ^ \X\, 
where TUW = {1,2,..., n} [161] . 

The LP bound for T-designs is formulated in a totally analogous manner, so we 
omit the details. This method provides lower bounds on the size of designs; see [52]. 
Due to the divisibility conditions on \ Y\ inherent in the definition of a T-dcsign Y in 
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the familiar cases, these bounds are not often as sharp as the corresponding bounds 
for codes, but in many cases, these are still the best known non-trivial bounds 

From now on, suppose that (X, TV) is the Hamming scheme H(n,q). The most 
traditional case here is that T is of the form {1, 2, . . . , 5— 1} for some 1 ^ S ^ n, so 
that we seek an upper bound on A q (n, 5), the maximum size of a code in X with 
minimum distance (at least) 5. Since Qij = Kj{i) is a (Krawtchouk) polynomial 
of degree j in i, the dual program can be stated entirely in terms of "Krawtchouk 
expansions" of polynomials: any polynomial / = ^2™=objKj satisfying (i) bo = 1; 
(ii) bj > (1 ^ j < n); and (hi) f(i) < (S < i ^ n), yields an upper bound 
A q (n,5) ^ /(0). Hence one may demonstrate feasible solutions for given ranges S 
without necessarily solving to optimality. Examples of bounds which can be derived 
in this way are the Plotkin bounds and the bound of McEliece et al. [124] . See 
also [59] for detailed discussions on the bounds of V. Levenshtcin. 

In 2001, A. Samorodnitsky [141] proved that, asymptotically, the optimum so- 
lution to Delsarte's LP bound is no better than the average of the upper bound of 
McEliece ct al. and the Gilbert-Varshamov lower bound: 

1 - H 2 (6) ^R^H 2 (^- y/6(l - efj (27) 

where H2(x) = — a;log 2 (a;) — (1 — x) log 2 (l — x) is the binary entropy function and 
R = limsup^^^ rT x log 2 A<i{n, On) is the asymptotic rate of the largest binary 
code with "relative minimum distance" 9 (= 5/n). In fact, these two bounds do 
not coincide for all < 8 < 1/2, so that even if the lower bound obtained from 
the Gilbert-Varshamov argument is close to the true optimal value of R, the linear 
programming method, alone, will never be sufficient to prove this. 

On the other hand, new ideas for obtaining upper bounds on codes are on the 
horizon. The theorem of A. Schrijver [143) applies semidchnitc programming to 
optimize over the positivc-scmidehnitc cone of the Tcrwilligcr algebra of H (n, 2); see 
i fTOl There have also been attempts to add more constraints to the program defining 
^lp using geometric arguments; see e.g., |127| §IV]. See also [127| I128| for another 
approach which focuses on the "holes" of codes in metric association schemes. It 
should be remarked that the recent determination of the kissing numbei0 in four 
dimensions (i.e., fc(4) = 24) [1321 1133] is also based on an improvement of the LP 
bound for spherical codes. See [17] . 

8. The ubiquity of codes and designs 

In this section, we summarize broad applications of Dclsartc Theory and the LP 
bound to various combinatorial objects. For the most part, these examples involve: 
an application; an association scheme (X, 1Z); a partially ordered set related to X; 
a partial order on the cigcnspaccs of (X, TZ) , somehow related to the first partial 
order. The short regular semilattices of Examples l5.1l and l5.2l are the paradigmatic 
examples, but we aim to emphasize more exotic settings of this same general flavor. 



^ 2 See, e.g., Table 4.44 in |93j where parameter sets for t-(v,k,X) block designs are ruled out 
(actually by Hacmers, Weug and Dclsartc) using linear programming. 

53 For S > (1 - l/q)n, the Plotkin bound uses / = K + b 1 K 1 where bi = l/(qS - n(g - 1)) to 
yield A q (n, S) ^ q&/{q$ — n(q - 1)). 

^The kissing number k(n) is the maximum number of unit spheres which can simultaneously 
touch the unit sphere in n-dimensional Euclidean space without pairwise overlapping. 
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In the cometric case, the eigenspaces are naturally ordered linearly; the Ei are in- 
dexed by the elements of the chain (totally ordered set) C n = {0, 1, . . . , n}. We may 
also consider products0 The eigenspaces of a product of m cometric association 
schemes are naturally ordered by a product of chains. Next, we consider exten- 
sion^3 of commetric association schemes. In this case, the partial order to consider 
on the eigenspaces is the quotient of the m-fold product of chains C n under the 
action of S m . This poset is obviously isomorphic to a downscl^ in Young's lattice 
(all partitions, or Ferrers diagrams, ordered by inclusion). It turns out that the 
theory allows us to take further products and extensions of all of these examples to 
obtain more. 

Since our main goal in this section is to demonstrate the widespread applicabil- 
ity of the theory discussed here, we now give a number of design-theoretic settings 
where Delsarte's concept of T-design applies. For each, we describe the combi- 
natorial objects in question, the association schemes in which they can be found, 
the relevant partial order on the eigenspaces of these schemes, and - relative to 
this indexing of eigenspaces - the subset T for which these objects are Delsartc 
T-designs. Applications are discussed in the references. 

Example 8.1. For block designs of strength t, i.e., t-(v, n, A) designs, the association 
scheme is the Johnson scheme J(v, n) with poset C n and T = {1, 2, . . . , t}. 

Example 8.2. For an orthogonal array of strength t, the ambient association scheme 
is the Hamming scheme H (n, q) with poset C n and T = {1, 2, . . . ,t}. 

Example 8.3. The incidence graph of a symmetric design is always distance-regular 
and has two Q-polynomial ordcrings. In |114j . a number of geometric substructures 
in finite projective spaces are shown to be Delsarte T-dcsigns in the corresponding 
cometric schemes. In all cases, the poset is C3 and we have T C {1,2}. 

Example 8.4. In [101j . Levenshtein studies systems of resilient functions with an eye 
toward cryptographic applications. The underlying combinatorial objects are split 
orthogonal arrays, which are Delsarte T-designs in a product H(n\, q) H(ri2, q). 
Here, the eigenspaces are indexed by C ni x C„ 2 and T = C tl x C t2 — {(0,0)}. 
Levenshtein was the first to derive the LP bound for these objects. 

Example 8.5. For some time, statisticians have been using mixed-level orthogonal 
arrays for experimental design, without knowing whether their constructions were 
as efficient as they could be. If such an experiment has factors with q\ possible 
levels (1 ^ i ^ m), then, in order to evenly test all i-tuples of factors, one seeks 
a Delsartc T-dcsign in the product ®™ =1 H(ni,qi). The poset on eigenspaces is 
x iL \Cn, and T = {{ji, . . . ,j m ) : < Y^iLiJi ^ The LP bound for mixed-level 
orthogonal arrays was derived by Martin [1 1 2j and, independently and simultane- 
ously, Sloane and Stufken [144] . 



^The product of two association schemes (X,7Z), (Y, S) with Bose-Mesner algebras A, B has 
vertex set X X Y and Bosc-Mcsncr algebra A ® B C C x xX ® C Y xY = C< x * y ) x xY), 

^^The Bose— Mesner algebra of the m-fold extension 1521 of an association scheme (X, 1Z) is the 
m tensor space of that of (X, TZ). The eigenmatrices of extensions are described in 11601 using 
generating functions and in 126 as Aomoto-Gclfand multivariate hypergeometric functions. See 
also [2iri25]. 

^In a partially ordered set (P, =^), a subset S C V is a downset (or lower ideal) if x S 5 and 
y =4 x always imply y <E S. 
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Example 8.6. One may consider the same set T for a product of Johnson schemes. 
An interesting special case is that of bipartite block designs (or mixed block designs) 
[111] , where points are colored with two colors, say red and white, and each block 
contains k\ red points and k 2 white points. For i + j ^ t, we require a constant 
number Ajj of blocks containing any i chosen red points and any j chosen white 
ones. If there arc v\ red and v 2 white points in all, then this is a T-design in 
J(vx,ki) g> J{v 2 ,k 2 ) where T = :0<i + j^t}. 

Example 8.7. A Room square of side n is a placement of all the unordered pairs 
of elements from = {1, 2, . . . , n + 1} into ("J 1 ) of the cells of an n x n array in 
such a way that each symbol appears once in each row and once in each column. A 
Room d-cube of side n is a d-dimensional cube of side n in which each 2-dimcnsional 
projection is a Room square of side n [51]. In the product scheme J(n + 1,2) <£> 
H(d,n), a Room <i-cube of side n is equivalent to a Delsarte T-design [112] of 
(minimal) size ( n J ) which is also an S-code where T = {(1, 0), (2, 0), (0, 1), (1, 1)} 
and S = {(1, 1), (2, 1), . . . , (1, d — 2), (2, d — 2)}. No example is known in which 
d> (n- l)/2. 

Example 8.8. Another recent application of Delsarte's theory of designs, and of the 
LP bound in particular, is the discovery of the ordered Hamming schem^^ [121j . 
The most important designs here are the ordered orthogonal arrays (OOAs), which 
in many cases give rise to (t,m, s)-nets |137[ II IT] which in turn provide quasi- 
Montc Carlo methods for numerical integration, optimization, and simulation. For 
an alphabet Q of size q, form relations Rq, R\, . . . , Re on Q by putting (x, y) € 
Ri if their maximal common prefix has length £ — ijfl The ordered Hamming 
scheme H (s, £, q) is the s-fold extension of the ^-class symmetric association scheme 
that results from this construction. Its eigenspaces are ordered by the downsct 
of all Ferrers diagrams in Young's lattice that fit inside a rectangle with s rows 
and £ columns. OOAs of strength t are characterized as Delsarte T-designs in 
H(s,£,q) where T contains all non-empty Ferrers diagrams with t or fewer cells. 
Since H (s, £, q) is self-dual, it was natural in [121j also to define "ordered codes" ; 
these turned out to be equivalent to the "codes for the m-mctric" studied in [139] . 
See also [1221I24]. 

We remark that Camion [37] also uses extensions of certain commutative asso- 
ciation schemes to derive Mac Williams identities for various weight enumerators of 
linear codes, such as complete (or spectral), Lee and split weight enumerators, all 
in a unified manner, together with the results of Sj6j 

Example 8.9. In the association scheme of the symmetric group S n , a most inter- 
esting class of T-dcsigns are the X-transitive sets of permutations [120] . Here the 
eigenspaces are indexed by partitions of n and the partial order is reverse domi- 
nance order <i. A set of permutations is A-transitive if it is a Delsarte T-design in 
this association scheme, where T = {fi : fi < A} — {(n)}. This association scheme 
also provides a framework for permutation codes (or arrays), and the corresponding 



Some authors refer to this as "NRT space," after Niederreiter, Rosenbloom and Tsfasman, 
whose earlier papers contained some properties of this association scheme without using the as- 
sociation scheme terminology at all. 

^In other words, we consider the £-fold wreath product H{1, q)l ■ ■ - l H(l, q) (see e.g., 11351 for 
a definition). 
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LP bound was studied in detail in [161j . These codes are of recent interest because 
of their application to data transmission over electric power lines; see |45j . 

The concept of "design systems" |115j is a far-reaching generalization of the 
poset structures associated with classical metric & cometric association schemes 
observed in e.g., [53] 1148) . and establishs a framework which includes all of the 
above types of Dclsarte T-designs as well as the most general bounds for themj^ 
It should be mentioned that there is another new approach to T-designs based on 
"coset geometries" [84] . This approach has the advantage that we can totally forget 
poset structures (though assuming transitive group actions), so that it may have 
the possibility to enable more flexible applications. 

9. The Terwilliger algebra 

Suppose that (X, TV) is a commutative association scheme. Fix a "base vertex" 
x G X. For each i (0 < i < n) let E* = E*{x), A* = A*(x) be the diagonal 
matrices in C XxX with (y, ^-entries (E*) yy = (A l ) xy , {A*) yy = \X\(Ei) xy . Note 
that E*E* = S lj E*, Yn=o E i = 7 > and moreover 

71 71 

xx x (28) 

fe=o i=o 

The E* and the A* form two bases for the dual Bose-Mesner algebra A* = A* (x) 
with respect to x. The Terwilliger (or subconstituent) algebra T = T(x) of (X,1Z) 
with respect to x is the subalgebra of C XxX generated by A and A* 163] 1164] I165j . 
The following are relations in T: 

E*A 3 E* k = iff p% = 0; E. t A*E k = iff g% = 0. (29) 

(The latter follows by computing the squared norm of EiA*Ek-) With the notation 
of Sj3l we also remark that if (X, TV} = G/K where K denotes the stabilizer of x in 
G, then T is a subalgebra of the centralizer algebra of tt\k- 

T C {M G C XxX : iv{g)M = Mir(g) for all g G A'}. (30) 

Equality in ([3H]) is known to hold for H(n, q) = (S q I S„)/(S q -i I S n ), for example; 
see [70]. 

Since T is closed under conjugate-transpose, it is scmisimplc and any two non- 
isomorphic irreducible T-modulcs in C x arc orthogonal. Describing the irreducible 
T-modulcs is an active area of research; see e.g. [1701 11081 [57] and the references 
therein. By (|2"9"|) we obtain 

Lemma 9.1 ([163]). Let W be an irreducible T -module. Then the following hold: 

(i) If (X,1Z) is metric with P-polynomial ordering {Ai}™ =0 , then A\E*W C 
E*_ t W + E*W + E* +1 W (0^i^ n), where E*_ x = E* +1 = 0. 

(ii) If(X,7V) is cometric with Q-polynomial ordering {-E;}™ =0 > then A\E^W C 
Ei_{W + E,W + E l+1 W (0<K n), where E_ x = E n+1 = 0. 



eo A 

design system attaches a poset <) to the eigenspace indices and embeds the vertex 
set X in a poset ("P, =<!) with an order-preserving surjection tp : (V, — » (J, <) satisfying three 
conditions. If M is the incidence matrix of X versus V (using =<! where XCp), then we require 
M to have no repeated columns and, for each j G J , the submatrix Mj obtained by restricting 
to columns x with <fi(x) = j must have constant column sum and column space Wj satisfying 
Vj C Wj C (Bi<jVi where Vi is the j th eigenspace of the scheme. 
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An irreducible T-module W is called thin (resp. dual thin) if dim E*W ^ 1 
(resp. dim EiW ^ 1) for all i. We remark that J(v, n) and H(n, q) are both thinW^ 
and dual thin, i.e., every irreducible T(a;)-modulc is thin and dual thin for every 
x G X . There are several infinite families of metric & cometric association schemes 
which have nonthin irreducible T-modules |1651 Note 6.2] such as classical forms 
schemes (e.g., bilinear forms schemes). In such cases, in general, the determination 
of all irreducible T-modules is yet to be settled, with the notable exception of the 
Doob schemes |153j . See also [3]. The irreducible T-modules of metric & cometric 
association schemes are often studied using the theory of tridiagonal pairs |85| . 
these being a generalization of Leonard pairs (Q. Namely, if (X, TV) is both metric 
and cometric then in view of Lemma 19.11 A\ and A* act on each irreducible T- 
module If as a tridiagonal pair (over C). We remark that W is thin (and dual 
thin) if and only if this tridiagonal pair is a Leonard pair. As of this writing, the 
classification has been worked out by Ito and Terwilliger [55J for the tridiagonal 
pairs over algebraically closed fields which have the most general "q-Racah" type. 
Their proof involves the representation theory of the quantum affine algebra Uq(sl%). 
See also [88]. It was earlier shown [87] that for the forms schemes there are four 
natural algebra homomorphisms from U q {s\2) to T, and that T is generated by 
each of their images together with the center Z(T)|fj It is also an important and 
urgent next step to "pull back" the above representation-theoretic information to 
the classification problem of metric & cometric association schemes. 

It would be a reasonable project to apply the progress on the Terwilliger algebra 
and the tridiagonal pairs to the analysis of codes and designs. This approach is 
still in its infancy, but it turns out that we may obtain several interesting results 
even from the elementary facts about T. As an example, we discuss the Assmus- 
Mattson Theorem [4], which gives a criterion as to when the supports of the words 
of a fixed weight k in a linear code in H{n, q) form a t-design (in J(n, k) U ~| 



Let W be an irreducible T-module. We recall the following parameters^ r = 
min{i : E*W ^ 0} (the endpoint); r* = min{z : EiW ^ 0} (the dual endpoint). Set 
1 = ^2 v eX y~ Then Ax = A*\, which is called the primary T-module. It is thin, 
dual thin and is the unique irreducible T-module in C x satisfying r — or r* =0. 

Suppose (X, 1Z) is cometric with Q-polynomial ordering {Ei}f =0 . A vector \ £ 
is a relative t-design with respect to x if Ei\ S CEiX for 1 ^ i $J t [54j . If 
(X, 1Z) is induced on the top fiber of a short regular semilattice (V,=4), then x is a 
relative t-design with respect to x if and only if for each u € V with rank(it) = t, 
J2 y ex u^yix^y} depends only on rank(ai A u) [5U Theorem 9.8]. In |158j . T was 
used to give a new proof of the following analogue of the Assmus-Mattson Theorem: 

Theorem 9.2 ( |54[ Theorem 8.4]). Suppose that (X,1Z) is cometric with Q-polyno- 
mial ordering {i?i}" =0 . Let Y C X be a code with characteristic vector \ and dual 



^This concept is in no way related to the "thin association schemes" of Zieschang 11771 . 

^ 2 These homomorphisms arise from an action of the q-tetrahedron algebra \El q I86| on T, so 
that their images are actually equal. 

^See 1821 for detailed discussions on the interaction of error-correcting codes with combinato- 
rial designs. 

^ 4 See footnote 1371 Though we do not (explicitly) use in this paper, but the following are also 
fundamental in the theory: d = \{i : E*W y£ 0}\ — 1 (the diameter); d* = \{i : EiW ^ 0}| — 1 
(the dual diameter). 
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distance S* . Set s x = \{i ^ : E*x 7^ 0}|. Then E* t x * s a relative (5* — s x )-design 
with respect to x for ^5 I n. 

Proof. Let U = (Ax)- 1 . Note that U is the linear span of all irreducible T-modules 
in C x with dual endpoint r* > 0. Set S = {i ^ : E*\ ^ 0}. Then 

where x|c/ denotes the orthogonal projection of x to U. Since A\ generates A* and 
takes s x (= \S\) distinct eigenvalues on J2jes ^jU, it follows that A*x|y is spanned 
by x\u,Alx\ii: ■ ■ ■ : (Al) s *'~ 1 x\u- Hence by Lemma RTT1 (ii) we find 

n 

A* X \u C J2 E * U - 

i=5* — s x + l 

This shows £"jA*x C CEii for 1 ^ z ^ 6* — s x , and the proof is complete. □ 

If the irreducible T-modules with dual endpoint at most S* — s x are dual thin 
in Theorem 19.21 then the conclusion can in fact be significantly strengthened: M\ 
is a relative (S* — s x )- design with respect to x for any M £ T . Note also that by 
dualizing the above arguments we may get another variant of the Assmus-Mattson 
Theorem for codes in metric association schemes. See jl 58} for the details. Theorem 
19.21 (as well as its dual) does not exactly coincide with the original when applied to 
H (n, q) with q > 2. It is interesting, however, to note that if (X, TV) is both metric 
and cometric then recent results on the displacement and split decompositions [170] 
can be successfully used to generalize the original version! 65 ! 

Theorem 9.3 ( |158j ). Suppose that (X,TV) is metric with P-polynomial ordering 
{^4i}™=0 an d cometric with Q-polynomial ordering {-E!j}™_ . Let Y C X be a code 
with characteristic vector x and dual distance S* . Set S x = min{i ^ : E*x 7^ 0}. 
Suppose t € {1, 2, . . . , n} is such that for every 1 $J r $J t we have 

\{r < % < n - r : EiX ^ 0}| ^ 6 X - r, or \{r < i < n - r : E*x ^ 0}| < 5* -r. 

If (X,7V) is induced on the top fiber of a short regular semilattice [V, =4), then for 
each M € T, X^gx u<y (-^X> V) * s independent of u =4 x with rank(u) = t. 

If (X,7Z) = H(n,q) and x is the zero vector (0,0, . . . , 0) (where € Q), then 
Theorem 19.31 shows that (the complements of) the supports of the words of fixed 
weight k in Y form a ^-design (in J(n, k) = J(n, n — k)) for every k. In particular, 
the conclusion of the original Assmus-Mattson Theorem is also true for nonlinear 
codes as 

A similar approach was also used in |158j to give a new proof of the minimum 
distance bound [113] for codes in H(n,q). We saw in ^6] that the MacWilliams 
identities for the weight enumerator of a linear code can be understood from the 
duality of A. The MacWilliams identities for the biweight enumerator |109j of a 

^Thc assumption on the semilattice structure is only for the sake of simplicity; sec 158 
Example 5.4]. It is assumed in 11581 Theorem 5.2] that the irreducible T-modules with endpoint 
at most t and displacement I170| zero are thin, but it follows from the results of 11511 that this 
condition is always satisfied. 

66 For example, the [12,6,6] extended ternary Golay code has covering radius three, and it 
follows from Theorem 19. 31 that a coset of weight three support 1-designs. 
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binary linear code can then be proved in terms of T for H(n, 2); see |118j . The 
harmonic weight enumerators of linear codes in H (n, q) and their Mac Williams 
identities studied in jHUHj use the harmonic analysis for the group S q -ilS n developed 
in J62[ [55] , so that we may view these as closely related to the theme discussed in 
this section. See also [34] [5] I154| for other proofs of the Assmus— Mattson Theorem 
based on harmonic analysis. 

We remark that the width and dual width of a code mentioned in Sj5] is quite 
compatible with the Terwilliger algebra theory. For instance, there is a more general 
approach [151| to the width, based on the Terwilliger algebra with respect to a code 
in metric association schemes. It is a generalization of the results of [7211167] on thin 
irreducible T-modules with e ndp oint one, and the width of a code and the tightness 
[92] of distance-regular graph^ are discussed together in the unified context of tight 
vectors; see also [80]. See |170| §8] for a generalization of Inequality (fT8]) . 

10. The semidefinite programming bound 

Throughout this section, suppose that (X, 1Z) is the binary Hamming scheme 
H(n, 2) = (S 2 I S n )/S n , so that X = Q n where Q = {0, 1}. Let x = (0, 0, ... , 0) 
be the zero vector and write T = T(x), E* = E*(x) (0 ^ i < n). Recall that T 
coincides with the centralizer algebra of K = S n acting on X. 

Let Y C X be a code. We consider two subsets IT , II 2 of G = S 2 I S n defined by 
IIi = {9 E G : x G gY}, Tl 2 = {g e G : x <£ gY}. For i G {1, 2}, let 

1 1 ' gen, 

where \gY G <C X denotes the (column) characteristic vector of gY. Since IT, II2 are 
unions of right cosets of G by K, it follows that Mg DP , M| DP G T. Moreover, since 
the XgY(XgY) T are nonnegative and positive semidefinite, so are Mg DP ,M| D p. By 
computing the inner products with the 01-matrices E*AjE£, we readily obtain 

-^sdp = ^v* E i A i E ki M snv = ^(^on - hjk)E*AjE%., 

i,j,k ijiifc 

where 

= \X\ \{(y,y',y")eY 3 :(y,y',y") satisfies (*)}| 
ljk \Y\ \{(y, y', y") G X* : (y, y', y") satisfies (*)}| ' 
and condition (*) is defined by 

(y,y')eRi, {y',y")&R h (y",y)&R k . (*) 

By viewing the Xijk as variables we get the following semidefinite programming (or 
SDP) bound established by A. Schrijver: 

Theorem 10.1 ([143]). Set 

4dp = ^sdp("-, S) = max^ f . J X 0ll 

subject to (i) A 00 o = 1; (H) < \jh < (Hi) \jk = K'j'k' if(i',f,k') is a per- 

mutation of(i,j,k); (iv) Ei^k^ikEZAjEl > 0; (v) Eij^ojj - h 3 k)E* A 3 E* k , 



°'Tight distance-regular graphs have many interesting combinatorial and geometric proper- 
ties, one of which is that every local subgraph is strongly regular with certain special nontrivial 
eigenvalues; see 1921 1721 . 
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0; (vi) Xijk = if {i,j, k} n {1, 2, . . . , S — 1} ^ (where > means positive semidefi- 
nite). Then ^(n, 6) ^ ^sdp- 

It is known that semidcfinite programs can be approximated in polynomial time 
within any specified accuracy by interior-point methods; see |172) . See also |69l 
§7.2] for a discussion on how to ensure that computational solutions do give valid 
upper bounds on A2(n, 5). While Dclsarte's LP bound is a close variant of Lovasz's 
$-bound, Schrijver's SDP bound can be viewed as a variant of an extension of 
the ^-bound based on "matrix cuts" |106j : see also [69l Chapter 6]. In fact, if 
we define a = (ao, Oi, . . . , a n ) by otj = Xouyf) (0 ^ i ^ n), then the condition 
that aQ is nonncgative is equivalent to the positive semidcfinitcncss of the matrix 
Mlp = J27=o ^ouAi, but since -Mlp = Mg DP + M| DP this is in turn a consequence 
of the positive scmidefiniteness of Mg DP and M| DP . A hierarchy of upper bounds 
based on semidcfinite programming was later proposed in |98j : 

i ( + > > ■ ■ ■ > *+° > ■ ■ ■ > Mn, 5). 

It turns out that ^lp = ^+ ^ ^sdp 5^ (■+ ' ■ Each of the £^ can be computed in time 
polynomial in n, but the program defining £^ already contains 0(n 7 ) variables. 
Two strengthenings of £sdp with the same complexity are also given in [98j . 

The SDP bound was also applied to the problem of finding the stability number 
of the graph (X, R n /2) for even n (known as the orthogonality graph) in [94], where 
it is shown (among other results) that for n = 16 the SDP bound gives the exact 
value 2304, whereas the LP bound only gives much weaker upper bound 4096. This 
problem arises in connection with quantum information theory [67]; see also |75j . 

As Mg DP , M| DP are 2™ x 2™ matrices, it is in fact absolutely necessary to simplify 
the program by explicitly describing the Wedderburn decomposition of the semisim- 
ple algebra T. The decomposition of T (as a centralizer algebra) was worked out 
in [62j in the study of addition theorems for Krawtchouk polynomials, but our 
discussion below emphasizes the use of T, based on [71] . 

Let W C C x be an irreducible T-module with endpoint r. Then W has dual 
endpoint r, and there is a basis {wi}™T r r for W such that 

Wi G E*W, A\Wi — (i — r + l)w.;+i + (n — r — i + l)w;_i (r ^ i ^ n — r) 

where u> r _i = w n — r +x = 0. Thus, the isomorphism class of W is determined by r. 
Moreover, it follows that 

(n — 2r\ 
■ _ r ) I \ w r\ I 2 {r <n- r). 

See [7T] for the details. The actions of the At on W may be described from the above 
information as the Ai are Krawtchouk polynomials in A\, but our argument goes 
as follows. For integers i, k, t such that ^ k ^ i ^ n and ^ t ^ min{fc, n — i}, 
we recall the following normalization of the dual Hahn polynomials found in [33] : 

Q^w»-(»i,)("r')'«(-*7T-r 1 1 

where \ k (z) = k(n — k) — z(n + 1 — z). If i + j + k is odd then E*AjE% = since 
H(n, 2) is bipartite, so suppose that i + j + k is even. Then it follows that 

E*A,E* k A 2 El = f3^ 2 E*A 3+2 E^+a^E;A J El+^ l : 2 E*A^ 2 El 
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where /3j' k 2 = (t + l)(i + 1 - k + t), Oj k = (k - t)(i - k + t) + t(n - i - t) and 
7^2 = (k+l-t)(n+l-i-t),witht= (j + k-i)/2. Using 2A 2 = A^-ri/ we find 
E^A^Wk = \ k (r)wk (r ^ k ^ n — r). Combining these facts with the three-term 
recurrence relation for the Q\' k |33|, Theorem 3.1], we obtain 

E*A jWk = Ql> k (\ k (r))E;A t _ k w k = (A fe (r)) Q ~ [ 

for r^k^i^n — r, ^ j ^ n such that i + j + k is even, where t = (j + k — i)/2. 
(The Qf' fe for i > min{fc, n — i} are formally defined by the recurrence relation |33[ 
Theorem 3.1].) Hence, after orthonormalization of the Wi, we get the following 
algebra isomorphism which preserves the positive-semidefinite cones: 

Ln/2J 

ip ■ T — > ^ j^(n-2r+l)x(n-2r+l) 
r=0 

where the r th block of f(Aj) is the symmetric matrix (a^'£)™^ r given by 
a?-' = a- = ( Afc W)(t» ("-D 1/2 (^)" 1/2 it i + j + k even, 



if i + j + k odd, 



for r < fc < i < n - r, < j < n. See also [1431 fi73]. 

The SDP bound has also been formulated for binary constant weight codes (i.e., 
codes in J(v,n)) in [143] and for nonbiliary codes in [701 EH]- The description of 
the irreducible T-modulcs becomes more complicated in this case, but this method 
turns out to improve the LP bound for many parameters. It seems to be an im- 
portant problem to decide whether it is possible or not to establish a suitable SDP 
bound for i-designs in J(v,n) or H(n,q). The SDP bound for spherical codes was 
formulated in [7J; it provides a new proof of fc(3) = 12 and fc(4) = 240 See also 

HZ). 
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